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Abstract 

The coherent contribution of all neutrons in neutrino nucleus scattering due to the 
neutral current offers a realistic prospect of detecting supernova neutrinos. As a 
matter of fact, for a typical supernova at 10 kpc, about 1000 events are expected 
using a spherical gaseous detector of radius 4 m and employing Xe gas at a pressure 
of 10 Atm. We propose a world wide network of several such simple, stable and low 
cost supernova detectors with a running time of a few centuries. 
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1 Introduction. 

The a typical supernova an energy of about 10° 3 ergs is released in 
the form of neutrinos [1],[2]. These neutrinos are emitted within 
an interval of about 10 s after the explosion and they travel to 
Earth undistorted, except that, on their way to Earth, they may 
oscillate into other flavors. The phenomenon of neutrino oscil¬ 
lations is by now established by the observation of atmospheric 
neutrino oscillations [3] interpreted as —» v T oscillations, as 
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well as v e disappearance in solar neutrinos [4]. These results have 
been recently confirmed by the KamLAND experiment [5], which 
exhibits evidence for reactor antineutrino disappearance. Thus 
for traditional detectors relying on the charged current interac¬ 
tions the precise event rate may depend critically on the specific 
properties of the supernova, in particular its distance from the 
Earth. This, of course, may turn into an advantage for the study 
of the neutrino properties. An additional problem is the fact that 
the charged current cross sections depend on the details of the 
structure of the nuclei involved. 

In recent years, however, it has become feasible to detect neutri¬ 
nos by measuring the recoiling nucleus employing gaseous detec¬ 
tors. Thus one is able to explore the advantages offered by the 
neutral current interaction. This way there are no problems as¬ 
sociated with uncertainties in the flux of any neutrino flavor due 
to oscillations. Furthermore this interaction, through its vector 
component, can lead to coherence, i.e. an additive contribution 
of all nucleons in the nucleus. Since the vector contribution of 
the protons is tiny, the coherence is mainly due to the neutrons 
of the nucleus. 

In this paper we will derive the amplitude for the differential 
neutrino nucleus coherent cross section. Then we will utilize the 
available information regarding the energy spectrum of supernova 
neutrinos and evaluate the expected number of events for all the 
noble gas targets. Then we will show that these results can be ex¬ 
ploited by a network of small and relatively cheap spherical TPC 
detectors placed in various parts of the world (for a description of 
the apparatus see our earlier work [6]). The operation of such de¬ 
vices as a network will minimize the background problems. There 
is no need to go underground, but one may have to go sufficiently 
deep underwater to balance the high pressure of the gas target. 
Other types of detectors have also been proposed [7],[8]. 

Large gaseous volumes are easily obtained by employing long 
drift technology (i.e TPC) that can provide massive targets by 
increasing the gas pressure. Combined with an adequate arnpli- 
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fying structure and low energy thresholds, a three-dimensional 
reconstruction of the recoiling particle, electron or nucleus, can 
be obtained. The use of new micropattern detectors and espe¬ 
cially the novel Micromegas [9] provide excellent spatial and time 
accuracy that is a precious tool for pattern recognition and back¬ 
ground rejection [10],[11]. The virtue of using such large gaseous 
volumes and the new high precision microstrip gaseous detectors 
has been recently discussed in a dedicated workshop [12] and 
their relevance for low energy neutrino physics and dark matter 
detection has been widely recognized. Such low-background low- 
energy threshold systems are actually successfully used in the 
CAST [13] solar axion experiment and are under development 
for several low energy neutrino or dark matter projects [6],[14]. 


2 Elastic Neutrino nucleon Scattering 


The cross section for elastic neutrino nucleon scattering has ex¬ 
tensively been studied. It has been shown that at low energies it 
can be simplified and be cast in the form: [1], [15]: 
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where is the nucleon mass and gy, g A are the weak coupling 
constants. Neglecting their dependence on the momentum trans¬ 
fer to the nucleon they take the form: 

1.27 

gv = < 2s\vl‘ 6 W + 1/2 ~ 0.04 , g A = —, (v,p) (2.2) 

1.27 

g v = -l/2,g A = —— , {y,n) (2.3) 


In the above expressions for the axial current the renormalza- 
tion in going from the quark to the nucleon level was taken into 
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account. For antineutrinos gA 


qa- To set the scale we write: 
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The nucleon energy depends on the neutrino energy and the scat¬ 
tering angle and is given by: 


T n = 


2 mN(E u cos 9)' 


{m N + E u ) 2 - ( E v cos 0) 2 
The last equation can be simplified as follows: 


T, 
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The above formula can be generalized to any target of mass m. 
It can be written in dimensionless form as follows: 


2 cos 2 6 
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The maximum energy occurs when 9 = 0 and depends on the 
neutrino energy (see Fig. 2.1). One can invert the above equation 



Wlrecoil 


Fig. 2.1. The maximum recoil energy as a function of the neutrino energy (both in 
units of the recoiling mass) 

and get the minimum neutrino energy associated with a given 
recoil energy. This is useful in obtaining the differential cross 
section (with respect to the recoil energy) after folding with the 
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neutrino spectrum. One finds: 


x = 


-1 + 


\ 


1 + 
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( 2 . 6 ) 


This is shown in Fig. 2.2 



"Irecooil 


Fig. 2.2. The minimum neutrino energy as a function of the recoil energy (both in 
units of the recoiling mass) 


3 Coherent neutrino nucleus scattering 


From the above expressions we see that the vector current contri¬ 
bution, which may lead to coherence, is negligible in the case of 
the protons. Thus the coherent contribution [16] may come from 
the neutrons and is expected to be proportional to the square 
of the neutron number. The neutrino-nucleus scattering can be 
obtained from the amplitude of the neutrino nucleon scattering 
under the following assumptions: 

• Employ the appropriate kinematics, i.e. those involving the 
elastically scattered nucleus. 

• Ignore effects of the nuclear form factor. Such effects, which 
are not expected to be very large, are currently under study 
and they will appear elsewhere. 
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• The effective neutrino-nucleon amplitude is obtained as above 
with the substitution 

q ^ j E N => 

with q the nucleon momentum and p the nuclear momentum. 

Under the above assumptions the neutrino-nucleus cross section 
takes the form: 



dT A J 

weak 


-[{iviy 


vl Hj v . 


+(Mv - M A )\ 1 




+ (Mi 


A-l T a 1 
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Where My and M A are the nuclear matrix elements associated 
with the vector and the axial current respectively and T A is the 
energy of the recoiling nucleus. The axial current contribution 
vanishes for 0 + => 0 + transitions. Anyway it is negligible in front 
of the coherent scattering due to neutrons. Thus the previous 
formula is reduced to: 



ZlMl (N 2 /4)F coh (A,T A ,E„), 

Z7T 

F coh (A, T a , E v ) = (l + + (l - 

(l _ A ~ 1 Ta 1 \ _ AtunTa 

V A m N E v /T a - l) E'l 


E„ J 

(3.8) 


The function F co h(A, T 4 , E v ) is shown in Fig 3.3 as a function of 
the recoil energy in the case of Ar and Xe for 10, 20,30 and 40 
MeV respectively. 


We see two reasons for enhancement of the cross section: 


• The overall A factor due to the kinematics, which is counter- 
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Fig. 3.3. The function F co h(A, Ta, E u ) as a function of the recoil energy Ta for, from 
left to right, E u = 10, 20, 30, 40 MeV. The results shown are for Xe on the left 
and Ar on the right 


acted by the smaller nuclear recoil energy when compared to 
the nucleon recoil energy for the same neutrino energy. This 
factor will be absorbed into the energy integrals, see the func¬ 
tion Ff 0 i d (A , T, ( T A )th ) below. 

• The N' 2 enhancement due to coherence. 


4 Supernova Neutrinos 


The number of neutrino events for a given detector depends on 
the neutrino spectrum and the distance of the source. We will 
consider a typical case of a source which is about 10 kpc, he. 
D = 3.1 x 10 22 cm with an energy output of 3 x 10 53 ergs with 
a duration of about 10 s. We will further assume that the en¬ 
ergy is shared equally by each neutrino flavor. Furthermore each 
neutrino flavor is characterized by a Fernri-Dirac like distribution 
times its characteristic cross section, i.e U v = 0.5 x 10 53 ergs per 
neutrino flavor, i.e. 


dN_ = El _ A x 4 

dE v u) l + exp{E u /T) JT 1 + e x 


(4.9) 


with J = A a constant and T the temperature of the emitted 
neutrino flavor. Each flavor is characterized by its own tenrpera- 
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ture as follows: 


T = 8 MeV for z/ /x , z/ r , z> M , F r and T = 5 (3.5) MeV for F e (y e ) 

The constant A is determined by the requirement that the dis¬ 
tribution yields the total energy of each neutrino species. 


U u = ^Jdx X 


U v 


J J"~ l + e x ^ A T 


Thus one finds: 

A = 0.89xlO 58 (z/ e ), 0.63 xlO 58 (z> e ) ,0.39xl0 58 (all other flavors) 


The emitted neutrino spectrum is shown in Fig. 4.4. 



E v MeV 


Fig. 4.4. The supernova neutrino spectrum. The short dash, long dash and contin¬ 
uous curve correspond to v e , u e and all other flavors respectively 


The differential event rate (with respect to the recoil energy) is 
proportional to the quantity: 

^W xFMA ’ Ta ’ xT) ^ (410) 

with A(T) = (0.89,0.63, 0.39) for z/ e , v e and all other flavors re¬ 
spectively. This is shown in Figs. 4.5 and 4.5. The total number 
of expected events for each neutrino species can be cast in the 
form: 




Fig. 4.5. The differential event rate as a function of the recoil energy T A , in arbitrary 
units, for Xe. On the left we show the results without quenching, while on the right 
the quenching factor is included. We notice that the effect of quenching is more 
prevalent at low energies. The notation for each neutrino species is the same as in 
Fig. 4.4 




MeV 


Fig. 4.6. The same as in Fig. 4.5 for the Ar target. 

No of events = C v {T)h(A,T, (T A ) th ), 

’ ,( A)th) ~ F fdd ( 40, T, 

with 




(4.11) 

(4.12) 


Ff 0 id(A,T, ( T A ) t h ) 


A 

1 


OO 


(T A )r 


C Ta ) 


dT A 

1 MeV 


x 


th 


, T a , xT ) 


x 


1 + e a 


(4.13) 


CAT) 


2n T 


9 



Where k is Boltzmann’s constant, P the pressure, V the volume, 
and Tq the temperature of the gas. 


Summing over all the neutrino species we can write: 


No of events = C v r(A) 


K(A, (T A ) th ) 
K{ 40, (T A ) th ) 


Qu(A) 


(4.15) 


with 


C v = 153 
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(4.16) 


In the above expression r(A) is a kinematical parameter depend¬ 
ing on the nuclear mass number, which is essentially unity. 

K(A, (T A )th) is the rate at a given threshold energy diveded by 
that at zero threshold. It depends on the threshold energy, the 
assumed quenching factor and the nuclear mass number. It is 
unity at ( T A )th ) = 0. The function r(A) is plotted in 4.7. It is 
seen that it can be well approximated by unity. 

From the above equation we find that, ignoring quenching, the 
following expected number of events: 

1.25, 31.6, 153, 614, 1880 for He, Ne, Ar, Kr and Xe(4.17) 


respectively. For other possible targets the rates can be found by 
the above formulas or interpolation. 

The function K(A, (T A ) t h) is plotted in Fig. 4.8 for threshold en¬ 
ergies up to 2keV. We see that the threshold effects are stronger 
in heavier systems since, on the average, the transfered energy is 
smaller. Thus for a threshold energy of 2 keV in the case of Xe 
the number of events is reduced by 30% compared to those at 
zero threshold. 

The quantity Qu{A ) is a factor less than one multiplying the to¬ 
tal rate, assuming a threshold energy {TjCjth = 100eV, due to the 
quenching. The idea of quenching is introduced, since, for low 
emery recoils, only a fraction of the total deposited energy goes 
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Fig. 4.7. The function r(A) versus the nuclear mass number. To a good approxima¬ 
tion r(A) ~ 1.0 (for definitions see text) 



Fig. 4.8. The function K(A, ( TA)th ) versus (TjCjth for various the nuclear mass num¬ 
bers without the quenching factor. From top to bottom He, Ne, Ar, Kr and Xe. (for 
definitions see text) 


into ionization. The ratio of the amount of ionization induced 
in the gas due to nuclear recoil to the amount of ionization in¬ 
duced by an electron of the same kinetic energy is referred to 
as a quenching factor Qf ac . This factor depends mainly on the 
detector material, the recoiling energy as well as the process con¬ 
sidered [17]. In our estimate of Qu(T A ) we assumed a quenching 
factor of the following empirical form motivated by the Lidhard 
theory [17]-[18]: 

Qfac(T A ) = n —777 J r i- 0.256 , r 2 ~ 0.153 (4.18) 

Then the parameter Qu(A ) takes the values: 

0.49, 0.38, 0.35, 0.31, 0.29 for He, Ne, Ar, Kr and Xe(4.19) 


11 



respectively. The effect of quenching is larger in the case of heavy 
targets, since, for a given neutrino energy, the energy of the re¬ 
coiling nucleus is smaller. Thus the number of expected events 
for Xe assuming a threshold energy of 100 eV is reduced to about 
560. 

The effect of quenching is exhibited in Fig 4.9 for the two inter¬ 
esting targets Ar and Xe. 



(T A )th —>MeV 


Fig. 4.9. The function K(A, (TjC)th) versus (TjCjth f° r the target Ar on the left and 
Xe on the right. The short and long dash correspond to no quenching and quenching 
factor respectively. One sees that the effect of quenching is less pronounced at higher 
thresholds. The differences appear small, since we present here only the ratio of the 
rates to that at zero threshold. The effect of quenching at some specific threshold 
energy is not shown here. For a threshold energy of 100 eV the rates are quenched 
by factors of 3 and 3.5 for Ar and Xe respectively (see Eq. (4.19). 

We should mention that it is of paramount importance to experi¬ 
mentally measure the quenching factor. The above estimates were 
based on the assumption of a pure gas. In our detection scheme 
the Xe gas carrier (A) is mixed with a small fraction of low ion¬ 
ization potential gas (B). Thus a part of the excitations produced 
on the Xe atoms could be transferred to ionization through the 
well known Penning effect as follows: 

A* + B—>A + B* + + e~ (4.20) 

Such an effect will lead to an increase in the quenching factor 
and needs be measured. 
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5 The NOSTOS detector network 


A description of the NOSTOS project and details of the spherical 
TPC detector are given in [6]. We have built a spherical proto¬ 
type 1.3 m in diameter which is described in [19]. The outer vessel 
is made of pure Cu (6 mm thick) allowing to sustain pressures 
up to 5 bar. The inner detector is just a small sphere, 10 mm 
in diameter, made of stainless steel as a proportional counter lo¬ 
cated at the center of curvature of the TPC. We intend to use 
as amplifying structure a spherical TPC [20] and developments 
are currently under way to build a spherical TPC detector using 
new technologies. First tests were performed by Filing the vol¬ 
ume with argon mixtures and are quite promising. High gains 
are easily obtained and the signal to noise is large enough for 
sub-keV threshold. The whole system looks stable and robust. 
The advantages of using the spherical detector concept are the 
following; 

(1) The natural radial focusing of the TPC allows to collect and 
amplify the deposited charges by a simple and robust de¬ 
tector using a single electronic channel to read out a large 
gaseous volume. The small size associated to small detector 
capacitance permits one to achieve very low electronic noise. 
In the present prototype the noise is as low as a few hun¬ 
dred electrons and has easily been obtained; with optimized 
low noise amplifiers we hope to lower it to the level of a few 
tenths. This is a key point for the obtaining a very low energy 
threshold, i.e. down to 100 eV, by operating the detector at 
moderate gain of about 100. Such low gains are easily ob¬ 
tained at atmospheric pressure and open the way to operate 
the TPC at high pressures. We target pressures as high as 10 
bar for Xenon gas. Even higher pressures by a factor 3-6 are 
aimed in the case of Argon gas in order to achieve, to first 
order, the same number of events. 

(2) The radial electric held, inversely proportional to the square 
of the radius, is a crucial point for measuring the depth of the 
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interaction by a simple analysis of the time structure of the 
detector signal. A position resolution of about 10 cm has been 
already obtained, a fact that is of paramount importance for 
improving the time resolution of the detector and rejecting 
background events by applying fiducial cuts. 

(3) Building a high pressure metallic sphere, for instance made 
out of stainless steel or copper, seems to assure an excellent 
quality of the gas mixture and turns out that a single gas 
filling with pure gas is sufficient to maintain the stability of 
the signal for several months. We are pushing the technology 
to improve the properties of the various elements in order to 
achieve stability over many years. 

(4) Big high pressure-secure tanks are under development by 
many international companies for hydrogen or oil storage, 
and therefore the main element of the TPC could be shipped 
at moderate cost. 

Our idea is then to build several such low cost and robust de¬ 
tectors and install them in several places over the world. First 
estimations show that the required background level is modest 
and therefore there is not need for deep underground laboratory. 
A mere 100 meter water equivalent coverage seems to be suffi¬ 
cient to reduce the cosmic muon flux at the required level (in the 
case of many such detectors in coincidence, a modest shield is 
sufficient). The maintenance of such system could be easily as¬ 
sured by Universities or even by secondary schools. Thanks to 
the simplicity of the system it could be operated by young stu¬ 
dents with a specific running program and simple maintenance 
every a few years. Notice that such detector scheme, measuring 
low energy nuclear recoils from neutrino nucleus elastic scatter¬ 
ing, do not determine the incident neutrino vector and therefore 
it is not possible this way to localize the Supernova. A cluster of 
such detectors in coincidence, however, could localize the star by 
a triangulation technique. 

A network of such detectors in coincidence with a sub-keV thresh¬ 
old could also be used o observe unexpected low energy events. 
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This low energy range has never been explored using massive de¬ 
tectors. A challenge of great importance will be the synchroniza¬ 
tion of such a detector cluster with the astronomical 7 -ray burst 
telescopes to establish whether low energy recoils are emitted in 
coincidence with the mysterious 7 bursts. 
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6 Conclusions 


In the present study it has been shown that it is quite simple 
to detect typical supernova neutrinos in our galaxy. The idea 
is to employ a small size spherical TPC detector filled with a 
high pressure noble gas. An enhancement of the neutral current 
component is achieved via the coherent effect of all neutrons in 
the target. Thus employing, e.g., Xe at 10 Atm, with a feasible 
threshold energy of about 100 eV in the detection the recoiling 
nuclei, one expects between 600 and 1900 events, depending on 
the quenching factor. We believe that networks of such dedicated 
detectors, made out of simple, robust and cheap technology, can 
be simply managed by an international scientific consortium and 
operated by students. This network comprises a system, which 
can be maintained for several decades (or even centuries). This 
is is a key point towards being able to observe few galactic su¬ 
pernova explosions. 
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